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Abstract 

A variational approach is presented to calculate the stress field gener- 
ated by a system of dislocations. It is shown that in the simplest case, 
when the material containing the dislocations obeys Hooke's law the vari- 
ational framework gives the same field equations as Kroner's theory. How- 
ever, the variational method proposed allows to study many other prob- 
lems like dislocation core regularisation, role of elastic anharmonicity and 
dislocation-solute atom interaction. The aim of the paper is to demon- 
strate that these problems can be handled in a systematic manner. 



1 Introduction 

Variational methods are frequently used in many fields of physics. They pro- 
vide a systematic and compact way to study different problems as well as allow 
relatively simple ways to incorporate various physical interactions. They make 
possible to introduce generalised coordinates and variables reflecting the sym- 
metry of the problem. Moreover, they are extremely important in numerical 
finite element simulations. In dislocation theory, however, variational methods 
are rarely applied. 

One of the few examples is the parametric dislocation dynamics elaborated 
by Ghoniem et at p,|. They have suggested a thermodynamics-based varia- 
tional method to establish the equations of motion for three-dimensional (3D) 
interacting dislocation loops. Their approach is appropriate for investigations 
of plastic deformation at the mesoscopic scale by discrete dislocation dynamics 
simulations. A fast sum technique for determination of elastic field variables of 
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dislocation ensembles is utilised to calculate forces acting on generalised coor- 
dinates of arbitrarily curved loop segments. Each dislocation segment is repre- 
sented by a parametric space curve of specified shape functions and associated 
degrees of freedom. 

Another example is the Phase Field method proposed by Rodney et al. to 
study dislocation dynamics They developed a general formalism for incor- 
porating dislocations in phase field methods based on the elastic equivalence 
between a dislocation loop and a platelet inclusion of specific stress-free strain. 
They also discuss how dislocations are elastically and dynamically coupled to 
any other field such as a concentration field. In their investigations special at- 
tention is paid to the treatment of dislocation cores after the discretisation of 
real and reciprocal space required by the computer implementation of any phase 
field method. In particular, they propose a method based on two length scales 
to account for dislocation cores much smaller than the grid spacing. 

In this paper we present a variational method that is quite different from 
the ones mentioned above. It is based on two assumptions: 

• The total deformation is the sum of the plastic and the elastic deforma- 
tions, i.e. we remain in small deformation limit. 

• The functional form of the Gibbs free energy-stress relation of the material 
under consideration is assumed to be known. 

First, the general form of the variational approach is outlined, then we discuss 
how dislocation core regularisation, elastic anharmonicity and dislocation-solute 
atom interaction can be treated within the variational framework. 



2 Derivation of the variational approach 

In order to study the collective properties of dislocations, as a first step, the 
stress field generated by a discrete dislocation system has to be determined. For 
a linear elastic medium the problem was formulated in a general way by Kroner 
[^. In this section we present a variational framework, which is equivalent 
with Kroner's theory for a medium obeying Hooke's law but it allows to study 
non-local and non-linear materials too. Elements of this variational approach 
have been presented in [H [S] , here we give a detailed derivation and clarify the 
meaning of the variational functional. 

It is commonly assumed that the final state of a body subject to shape 
changes is reached by two subsequent steps, a plastic and an elastic deformation 
[6], i.e. 

where Uj is the displacement vector, /3y , /3?. are the total, elastic and plastic 
distortions, respectively. In our analysis we remain in the small deformation 
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limit, i.e. we neglect the third term in the right-hand-side of the above equation. 
Thus the total deformation 



=(/5|,+/3,.)/2 (2) 



is the sum of the elastic eij — {(3ij + f3ji)/2 and plastic ef^ = {(3fj + 
deformations: 

The elastic deformation tij is generated by the stress field aij . Our aim now 
is to formulate the condition for the stress in the presence of a given plastic 
deformation field. For that purpose a convenient way to proceed is to write 
Cij as the functional derivative of the Gibbs free energy G[a\, functional of the 
stress tensor components Uij, i.e. 

With the above relation equation ^ reads as 

We shall use below the incompatibility operation —eikmSjhidkdiTmn on a tensor 
Tmn, where Cikm denotes the antisymmetric unit tensor (sometimes it is defined 
with opposite sign). Since the incompatibility of e*„„ is identically zero, by 
taking the incompatibility of equation ([5|) one obtains 

eikm,ejindkdi- = r]ij, (6) 

0<Jmn 

where 

Vij ■= eikmejindkdieP^n (7) 
is the incompatibility tensor of the plastic deformation. In terms of the dislo- 
cation density tensor , it has the form [3] 

Vij = 2 i^3irn.diaim + eumdiajm) ■ (8) 

So, we ended up with an equation, ([6]), which relates the stress field to the 
dislocation distribution and can be used to determine the former if the latter is 
given. 

We should add the condition of equilibrium, djaij = 0, which is automat- 
ically fulfilled by representing the stress tensor as the incompatibility of the 
stress potential Xij 

The functional derivative by the stress potential x is easily expressed via partial 
integration from 

5G = j dV-^-5(Jij ^ J dV-^^eikmejindkdiSxmn = J (10) 
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to yield 

eikmejindkdi- = - — , (11) 

damn OXij 

where G is considered as a functional of the tensor Xij only through its de- 
pendence on aij by the definition Finally, by the use of equation (|lip the 
equilibrium condition ([5]) becomes 



SG 

By introducing the 'plastic free energy' as 



(12) 



P[x,ri]--=G[x]+ UVl^^,X^J. (13) 



equation (|T2|) can be reformulated as 



0. (14) 



This way we arrived at a variational principle which yields the equilibrium 
equations for the stress potential in the presence of dislocations. The functional 
P has a simple interpretation in that G\x\ is the elastic Gibbs free energy in 
the absence of defects, complemented with the term representing the interaction 
between dislocations and the stress field. Without dislocations, the first term 
alone is the functional to be extremised, yielding the bulk equilibrium equations. 
In our formulation the incompatibility tensor assumes the role of the 'charge' of 
dislocations, that interacts with the stress potential linearly, and through the 
variation by the potential it enters the equilibrium condition. 

A further interesting relation can be obtained by taking the plastic free 
energy (fTS]) at the equilibrium stress potential x^'^- Substituting equation ^ 
into equation ()13|) and performing partial integrations one concludes that 

P[X^^ ^\ = G[a^^] J dV^-^^a^^l (15) 

which is just the Legendre transform of the Gibbs free energy and is thus the 
elastic free energy of the system as 

F[rj]^P[x"',7j]. (16) 

Note that, before extremisation, P[x, ">]] is, in general, not the free energy, but 
it takes the value of the free energy in elastic equilibrium, when it becomes a 
functional of the dislocation distribution through the incompatibility tensor r]. 

To summarise this section, we can state that we need two constitutive rela- 
tions to determine the stress field generated by a set of dislocations. The first 
one is given by equation ^ expressing that the final state of a body subject 
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to plastic deformation is reached by a plastic and a subsequent elastic defor- 
mation. In our considerations we neglect the cross term between plastic and 
elastic deformations, i.e. we remain in the small deformation limit. The other 
constitutive relation we need is the Gibbs free energy versus stress functional. 
With these the stress state generated by a dislocation system can be determined 
by finding the extremum of the functional given by equation (jl3p . It should be 
mentioned, that the approach explained above provides only the bulk equations. 
To get a unique solution, one has to set up boundary conditions appropriate for 
the problem considered which in most cases cannot be naturally obtained from 
the variational approach proposed. Since in this paper we concentrate only 
on bulk properties, in the problems explained below we considered either infi- 
nite size systems with vanishing fields at infinity or we used periodic boundary 
conditions. 

Formerly, we have used previous versions of the variational formalism to 
treat the effect of coarse graining and the correlations it entails, and recognised 
its use to derive dynamical equations for dislocation densities via the phase field 
method. Further studies in that direction will be presented in other publica- 
tions. In the following, we discuss how the extremum principle obtained can 
be applied in a variety of equilibrium problems, like dislocation core regularisa- 
tion, studying anharmonic effects and the interaction of dislocations with solute 
atoms. 



3 Local linear medium 

Wc first demonstrate the variational principle on the elementary example of a 
local, linear material obeying Hooke's law. In this case the Gibbs free energy is 
a quadratic functional of the stress as 

GoM:=-/i..,%.a.dF, (17) 

where Sijki is the elastic compliance tensor. Hence, the plastic free energy given 
by equation ([13]) reads as 



dV. (18) 



One can find from the equations (|T4)l and p8)) that the stress function fulfils a 
system of coupled, linear, fourth order partial differential equations 



eksteiuzdsdudodqSijkiXpr = Vtz- (19) 



Of course, the left-hand-side is just the incompatibility of the elastic deformation 
tensor. The above equation greatly simplifies for isotropic materials with shear 
modulus /i and Poisson's ratio v, when it becomes (see [5]) 

A'x:, (20) 
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where 

if x'ij satisfies the 'gauge condition' 

^^x'^, = 0. (22) 

3.1 Plane problems 

In the rest of the paper we restrict our analysis to plane (2D) problems, i.e. we 
consider only systems of straight dislocations extending parallel to the z direc- 
tion. After a long but straightforward calculation one can find that for edge 
dislocations the plastic free energy functional reads as 

1 



-(Ax) + x(52a3i ~ dia32) 



d'r, (23) 



4^l 

where x •= X33 is now a single component stress function (the other components 
of Xij vanish) and the stress tensor components are 

0-11 = -d2d2X, 0-22 = -didix, (J12 = did2X- (24) 

In the plastic free energy above the incompatibility tensor rjij has been expressed 
by the dislocation density tensor aij, which now has only two components. The 
extremum condition SP/Sx ~ leads to the fourth order partial differential 
equation as 

-^A^X = 52a3i - 9iq;32. (25) 

Although later in this paper we do not consider screw dislocations, we sum- 
marise their case as well. To our knowledge this was the first application of a 
variational principle in terms of the stress potential by Berdichevsky [7] , which 
motivated us to look for a variational description of general, including edge, 
dislocations. The plastic free energy now is 

P[0,a] =y (^-^|V0p + </,a33^dV (26) 

with 

<j) -.^-81X23 + 82X31, (27) 
and the relevant stress components are 

(T23 = ~di(f>, (Ti3 = 92 0. (28) 

The corresponding extremum condition gives the Poisson's equation 

-A(j) = -ass- (29) 

^^ 

It has to be mentioned that the above equations obtained by the variational 
approach are certainly equivalent with the ones derived earlier by different meth- 
ods (see Kroner [3]). In this section we just demonstrated how the variational 
method works for a classical local linear medium. 
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4 Dislocation core regularisation 



The significance of dislocation core regularisation is widely known. It is not only 
necessary for the explanation of physical core effects, but also to eliminate sin- 
gularities in a physically well founded manner in numerical simulations. There 
are many different propositions for dislocation core regularisation (for a review 
see [5]) but, as it is explained below, the variational approach offers a natural 
way to regularise the singular stress at the dislocation line. 

It is common in phase field theories that surface or size effects are captured by 
introducing appropriate 'gradient terms' in the energy functional. The concept 
can be applied in dislocation theory too, but as we have recognised above, the 
physical properties of a material are determined by the functional form of the 
Gibbs free energy-stress relation. So, the gradient terms should be introduced 
in the Gibbs free energy and we have to add terms that depend on the gradient 
of the stress. In first order linear approximation one can consider the 'non-local' 
Gibbs free energy 

G_.e.M:=GoH-6^/A^.,.™„(a.a,.)a..„^ (30) 

where Nijkimn is a constant tensor with inverse stress dimension and b is the 
Burgers vector. (6^ is separated from Nijkimn to indicate the relative order 
between Go and the gradient dependent term.) From Gnon-iocai [c] the corre- 
sponding PiXji]] ti^s to be constructed as it is explained in section [51 

It should be mentioned that non-locality could be introduced on a much 
more general way by taking the Gibbs free energy in the form 

Gnon-iocaiH ■= - J ^<T^j (r)5„ w (r " r')akiir')dVdV' (31) 

where Sijki{r) is a function which goes to zero fast enough if |r| — > oo, but in a 
first order approximation, if its range is of the order of the lattice constant, it 
obviously gives the same as (pO]) . 

To demonstrate how the non-local term introduced above results in disloca- 
tion core regularisation let us consider a single straight dislocation embedded in 
an isotropic medium. From equations ([23)1 and ([30)1 for a single edge dislocation 
at the origin 

P[X] = J [(^X)' + a'lVAxH + x92<5(r)| dV, (32) 

where a is a parameter with length dimension that is in the order of the lattice 
constant. Here, for the sake of simplicity, we considered only the simplest pos- 
sible isotropic gradient term from (|30p but the general case can be treated in a 
similar way. The corresponding equilibrium equation has the form 

A\~a^A^X^P^d2S{r). (33) 
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The above equation can be solved numerically quite easily. By taking its Fourier 
transform one can find that 



iqy 



from which, according to equation (j24p . the Fourier transform of the resolved 
shear stress reads as 

By numerical inverse Fourier transformation one obtains the shear stress plotted 
in figure [T] To demonstrate the difference more explicitly the shear stress along 
the X axis is plotted in figure[2] As it can be seen, for x larger than about 10 a, 
is close to the classical stress field ai2 oc 1/r. So, as it is expected, the 
'gradient term' introduced above influences only the central core region. 
In a similar way, for screw dislocations from equations ()26|) and (I30p 

P^J[-^ [|V0P + c2(A0)2] + 60^(r)| dV, (36) 

where c is a constant. Such approach was also made in [7]. One can obtain that 
the extremum condition 5P/ 5(j) = is fulfilled if satisfies the equation 

A(t>-c^/\'^(t> = -h^i&{r). (37) 

Like for edge dislocations, numerical solution of the above equation shows that 
the second term of the left hand side of the equation results in stress regulari- 
sation in the vicinity of the dislocation line. 

To summarise this section, we have shown that by adding appropriate terms, 
depending on the gradient of the stress tensor, to the Gibbs free energy func- 
tional one can regularise the dislocation core. We have discussed in detail only 
the simplest possible 'isotropic' core regularisation, but one can study more gen- 
eral cases too. Comparing the present approach with other methods suggested 
by Aifantis et al. jHUHji and Lazar [TD], although one can find some similarities, 
but the major difference is that in these works the core region is regularised by 
spreading out the dislocation density tensor atj , while in our analysis the dislo- 
cation density tensor remains proportional to a Dirac delta, like in the classical 
non-regularised case. 

Finally, it should be stressed, that the physical justification of the results 
obtained above requires further detailed investigations. This is beyond the scope 
of this paper, our aim was only to demonstrate that the problem can be treated 
within the variational approach in a straightforward manner. 



5 Role of elastic anharmonicity 

There is quite an extended literature on how to deal with dislocations in the ge- 
ometrically and materially non- linear setting [TTl [T^l [13] . Among other things 
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X [a] 

Figure 2: Shear stress along the x axis obtained with and without core regular- 
isation. (G := A*/[27r(l - v)\) 

the stress field of single straight edge and screw dislocations are analytically 
determined. So far, however, the interaction energy of dislocations is not cal- 
culated for the non-linear case. Recently, Nabarro and Brown have suggested 
that the energy difference between dipoles with interstitial and vacancy types 
can play an important role in the physical properties of PSBs [Ml [15]. As it was 
realised by them, in order to account for this energy difference one should go 
beyond linear elasticity. Allowing a certain quadratic term in the stress-strain 
relation they were numerically able to estimate the energy difference between 
the two types of dipoles [15]. In this section we show that the problem can be 
treated in a systematic manner within the variational approach explained in 
section 2. 

Before we start our analysis, we have to compare the different possible 
sources of non-linearities in dislocation theory. As it is discussed at the be- 
ginning of section 2, according to equation ([1]), the general form of the total 
distortion contains the product of elastic and plastic distortions. This term 
was neglected in the subsequent considerations. However, as it is explained be- 
low, elastic anharmonicity plays a much more important role. Up to quadratic 
terms anharmonicity can be accounted by introducing a quadratic term in the 
stress-strain relation: 

(38) 
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where Lijki and Kijkimn are the first and second order elastic moduH, respec- 
tively. According to the MD simulations of Chantasiriwan and Milstein [16] , for 
cubic metals the ratio of the first and second order moduli is (5 = i^im/Ln ^ 
— 10. {S = —9.8 was reported experimentally on steel by Sommer et al. [T7].) In 
contrast with this, in equation ([1]) the 'coupling constant' between the linear and 
the quadratic terms is 0.5. Due to this, the non-linear effects related to elastic 
anharmonicity are an order of magnitude higher than the ones connected to the 
other non-linearity mentioned above. So, the assumptions needed to derive the 
variational framework are justified. 

Since in the variational method proposed the Gibbs free energy-stress re- 
lation represents the physical properties of the material under consideration, 
anharmonicity has to be taken into account by adding cubic and/or higher or- 
der terms to the quadratic Go[a] relation given by equation pT)) : 

G'anharm 

[a] Go[(t] + J [Gij 

klmn ^ij ^kl ^mn 

+ 0(4)]dV^, (39) 

where Cijkimn is a constant tensor of rank six. In the following, we briefly show 
that the anharmonicity leads to an extra term in the interaction energy of two 
parallel edge dislocations. The reader can find the detailed derivation in [5]. 
For simplicity we consider isotropic medium and take the following Gibbs free 
energy 



G[x] ■■= J 



2 



-^(Axr+e^—iAx)' 



dr, (40) 



where D := (1 — and fa > is a dimensionless coupling parameter. 

Since, according to equation Ax = ~{(Jii + (J22) = 2p, the anharmonic 

term given above is proportional to the third power of the pressure p. So, it 
represents 'inclusion-exclusion' asymmetry. One can find that the other possible 
cubic terms in ([39]) result in quite similar dislocation-dislocation interactions 
than the one considered here. 

The corresponding plastic free energy defined by equation (|13|) is 



r, 



P[X]= J {^-j{Axf + e,^{Axf+x[bidySir-r^) + hdy6ir-r2)]jd^ 

(41) 

where bi = (61 , , 0) and 62 = (62 , , 0) are the Burgers vectors ( | &i | = 1 62 1 = : b), 
and ri and r2 are the positions of the two dislocations, considered as fixed for 
now. 

In order to calculate the interaction energy between two dislocations, we 
make use of the property that the plastic free energy extremised by the stress 
function will assume the value of the free energy. So, if we start out from the 
plastic free energy with two dislocations in fixed positions, extremise, then sub- 
stitute the equilibrium stress function, the term that depends on both positions 
will be just the sought interaction energy. While the leading pair energy from 
linear elasticity is long known, the anharmonic correction will be determined as 
follows. 
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First, we have to calculate the stress function x(r, ri, generated by two 
dislocations. One obtains from equations (fT4)) and ((4T|) that 



DA\ ~ e^D^A [(Ax)'] - bidy5{r - n) + hdySir ~ rs). (42) 

Since Ca is a small parameter, up to first order in ea, the solution of equation 
can be found as 

X{r,ri,r2) = Xi(r - ri) + X2(r - ra) + ea^(T-, ri, ra), (43) 

where 

Mr):=jXoir), i = l,2 (44) 
is the solution of the harmonic single dislocation problem 

DA\o - bdyS{r), (45) 

Xo(r) = ^Vln(r). (46) 

After substituting equations ([^S)) and (HS)) into equation one can find that 
if ea is small than ^ satisfies the equation 

A^ = D{Axi + Ax2f- (47) 

According to the general expression the plastic free energy (|¥T|) at the 
extremising stress function gives the free energy. Substituting equation (j42p 
into (PT|) gives 



F = P[x] = 1 



2 



-(Ax)'-ea — (Ax)3 



dV, (48) 



wherein we now understand the equilibrium x field. Up to terms linear in the 
coupling parameter ea, the total free energy is 

^^ = Fo + eai^i, (49) 
where Fq is the well-known energy in harmonic approximation and 

Fi:= y |i^A(xi+X2)Ae-^[A(xi+X2)]'|dV (50) 

From equations (gT]) and (|5n)) 

= J ^[A{xi + X2)?d\. (51) 

This can be rewritten as 

F^^D' J [(Axi)(Ax2)' + (Ax2)(Axi)']dV + C, (52) 
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where C is a constant representing the correction to the single dislocation energy. 
While this is singular in the present approximation, it does not contribute to the 
dislocation-dislocation interaction force. Thus, the Vi(r) :— Fi — C interaction 
energy correction has the form 



where r := ri — r2 is the relative coordinate of the two dislocations. 

The integral in the above energy correction can be calculated analytically. 
Since it is not straightforward, we give the detailed calculation in the Appendix. 
One finds that 



where Tq is the inner cut-off radius. Please notice that a term is missing from 
the energy correction form (and so the extra interaction force) given in jSJ. 
It is useful to summarise some of the properties of Vr. 

• If the Burgers vectors of the two dislocations are the same, Vi vanishes, 
so the non-linearity considered generates extra interaction between dislo- 
cations of parallel Burgers vectors only if their signs are opposite. 

• Vi decays as ln(r)/r with a complicated angular dependence. 

• Vi is antisymmetric in y. 

• Vi is symmetric in x. 

• If we take a dislocation and then we put another one with opposite sign 
above or below, the interaction energy of the two configurations arc not 
the same due to the extra interaction term Vi. So, as it is expected, the 
anharmonicity introduced recovers the well know fact that the interaction 
energy of the interstitial and vacancy type dislocation dipoles are not the 
same. 

A remarkable consequence of the extra term in the dislocation-dislocation 
interaction energy is that it strongly modifies the relaxed dislocation morphology 
developing from an initially random dislocation configuration. 

As it can be seen in figures [3] and [H in the harmonic case the relaxed configu- 
ration consists of randomly placed dipoles and short walls, while the presence of 
the extra interaction term leads to the formation of dense multipoles arranged 
in a structure with some periodic character [5]. 

6 Dislocation— solute atom interaction 

Solute atoms can strongly modify the collective properties of dislocations. Among 
other things they can lead to plastic instabilities (for a recent review see [T5]V 




(53) 




(54) 
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Figure 3: A relaxed dislocation configuration with Ca = (harmonic case). 
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Figure 4: A relaxed dislocation configuration with ea = 5 • 10 "^L/b, with L 
being the system size (anharmonic case). 
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In this section we show that the effect of solute atoms can be easily incorporated 
into the variational framework of section [2l We restrict our considerations for 
straight edge dislocations with Burgers vectors parallel to the x axis, but it is 
straightforward to generalise the method for 3D. 

The Gibbs free energy of a coupled system can be always given as the sum 
of the Gibbs free energy of the two uncoupled systems and a coupling term. 
Therefore, if we add to the plastic free energy the Gibbs free energy contribution 
of the solute atoms than, after extremisation in the stress function, we shall 
arrive at the free energy of the dislocation-solute system. According to equation 
the plastic free energy of the parallel edge dislocation system considered 

is 



dV, (55) 



4m 

where k is the signed dislocation density (geometrically necessary density, GND), 
defined as k := a^i/h, that is the only non- vanishing component of the disloca- 
tion density tensor for the dislocation geometry considered. 

For the solute atoms we assume that their concentration c is close to the 
equilibrium concentration Coo- In this case the Gibbs free energy of the solute 
atoms can be given with the quadratic form 

Gc[c] j a{c-c^fA\, (56) 

where a is a constant (which may depend on Cqq). 

To determine the form of the coupling term we use the well-known fact that 
beside concentration gradient, the pressure gradient also causes solute atom 
diffusion. According to the principles of irreversible thermodynamics, the solute 
atom current is proportional to the gradient of the chemical potential /i^ = 
bG 1 8c. So, taking the coupling term in the form 



G, 



couplinj 



;[x,c] := y /3cpd^r, (57) 



where /3 is constant, the total plastic free energy 

^'[X, C] = Pd[x, + Gc[c] + Gcoupling[x, c] 



i-^(Ax)' + 6x(92k) + a(c - Coo)' - /3cAx 
4/i 



d r 



(58) 



obviously gives the right free energy together with the expected form of [i^. 
It should be mentioned that ([55)1 results in linear equations for the stress and 
the solute atom concentration. Non-linearity can also be treated within the 
framework proposed, but it is out of the scope of the paper. 

Since this is mostly a demonstration of the way the variational approach 
works, we now restrict our analysis only to static problems. Dynamical as- 
pects of the interaction between dislocations and solute atoms will be consid- 
ered in forthcoming publications. One can find from the equilibrium conditions 



15 



5P/5x ^ and 5P/5c = that 



1-v 



A2 



(59) 



2/i 



and 



a(c - Coo) = /?Ax. 



(60) 



By combining the two equations we get that 



l-v 0- 
V — 

2/i Q 




(61) 



A remarkable feature of the above equation is that, apart from a constant mul- 
tipher, the functional form of the stress function x is not effected by the solute 
atoms. Moreover, the solute atom concentration is proportional to the pressure 
caused by the dislocations. Certainly, the result obtained is not new, it is the 
well-known Cottrell atmosphere of solute atoms around dislocation lines [12]. 
But it illustrates very well the fact that the coupled system of dislocations and 
solute atoms can be treated with the variational framework suggested in sec- 
tion [2l However, it is important to mention that the collective dynamics of the 
coupled system can also be studied with the variational method. This will be 
discussed in a forthcoming paper. 

7 Conclusions 

A variational framework has been presented to determine the stress field and the 
interaction energy of dislocation systems. It is demonstrated that the method 
allows to study core regularisation problems, the role of anharmonicity and 
dislocation-solute atom interaction. Furthermore, Groma et al. have discussed 
in detail how the dynamics of coarse grained dislocation densities can be ob- 
tained from an effective variational free energy [2DJ. This recognition led to 
the theory of Debye-like screening of dislocations in equilibrium [?] and during 
relaxation [5T]. Certainly, like in other fields of physics, one can obtain the 
same results without using the variational approach, but this offers a system- 
atic way for the treatment of a wide range of different problems. It should be 
also mentioned that other variational principles can also be constructed, which 
give the same result [22]. The advantage of the one presented here is that the 
variational functional is minimised with respect to the stress potential. So, it 
gives equations directly for the stress function from which the stress can be 
calculated directly. 
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Appendix. Calculation of the anharmonic poten- 
tial 



As it is explained in section 5, elastic anharmonicity leads to the extra dislocation- 
dislocation interaction energy (|53p : 



V^{r) 



-A 



y-y 



^ rd^r'. 



(62) 



where 



A := 



{bi - b2)b^ 



(63) 



The integral can be evaluated analytically leading to equation now we 
give the derivation. 

Let us first analyse the asymptotes of the integrand. At infinity it decays 
like and near r ^ diverges like l/|r' — r|, so in these cases integration 

gives finite result. Contrastively, near the origin the integrand is proportional 
to l/r'^, so a lower cut-off rg has to be introduced and the integral will exhibit 
a logarithmic singularity. This is an artificial singularity, however, because it 
comes from the core region, where interactions are regular in reality. The proper 
treatment would be by core-regularised anharmonic analysis, which is relegated 
to later work. 

To proceed with the calculation, we notice that ([5^ is the convolution of 
two functions, which can be rewritten as follows 



y_ 



1 
2 

dy 



d 



Qy2 ( ) j,2 

In (r) , 



(64) 
(65) 



where the lower cut-off tq is omitted from the logarithms but will be reinstated 
in the end. Substitution into ((62)) yields 



Vi{r)^^[V,^{r) + V,2{r)], 



with 



Vii{r) dy\n{r) ® dl\n{r) 



ln(r) (g)5^1n(r), 



Vi2{r) := -dy ln(r) (g) r ^ = - ln(r) dyV' 



where ® denotes convolution: 



(66) 



(67) 
(68) 



/(r)®5(r) := / A\' f{r - r')g{r'). 



(69) 



We then recall that the radial logarithm is the Green function of the Poisson's 
equation in the plane, so a convolution like 4){r) = ln(r) g{r) is just the 
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solution of A(/)(r) — 2TTg{r). Therefore, to calculate (|67p and ([55]) . we have to 
solve the Poisson's equation with the inhomogeneous terms dy ln(r) and —dyr~^. 
But since those functions are derivatives of radial functions, it suffices to solve 
the radial Poisson's equation and differentiate the solution later. The radial 
solutions needed for our considerations are 

g(r)-ln(r) => 0(r) = V [ln(r) - 1] , (70) 

g{r) = -r^'^ (/)(r) = -tt In^(r), (71) 

whence we get 

V^^ir) = ^^'ylnir), (72) 

Vi2{r) = -Trdy\ii^{r). (73) 

Reintroducing the cut-off tq into the logarithms and performing the differenti- 
ation gives by and ([55)1 the leading anharmonic term to the dislocation- 
dislocation interaction energy: 



Viir) = 



(&i - b2)b' 



y{3x^ +y^) r,y , f r 



H V ro 
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